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NOMENCLATURE 

a 

n 

- nth  complex  coefficient  In  Fourier  series  expansion  of  F(0) 

1 

* 

A 

- wave  amplitude 

A^iojT) 

- Fourier  transform  of  r|^(t) 

B [ © ; u>|  6q; 

ii)Q]  - beam  pattern 

C i C 

^ j - covariance  matrix  (spatial) , matrix  elements  ' 

e{  } 

- expected  value  operator  (statistical  average) 

E (0) 

m 

- equivalent  to  exp[j  k*r  ] 

— — m 

F(0) 

- normalized  directional  spectrum  spreading  function 

f 

- frequency  in  cycles /second 

F.O.M. 

- figure  of  merit 

H 

- wave  height  (twice  wave  amplitude) 

V > 

- kth  order  Bessel  function 

k,  k(w) 

- wavenumber  (vector),  wavenumber  (scalar) 

n(x,t) 

- noise  signal 

N 

- number  of  array  elements  or  sensors 

N(u>) 

- noise  power  level  in  frequency  domain 

p(0o;a,o) 

- weighted  measure  of  array  performance 

7(a,o) 

- average  of  P(6  ;u>  ) over  all  0 angles 
o o o 

Q' 

- figure  of  merit,  F.O.M. 

m 

Vtt) 

- element  of  the  inverse  of  covariance  matrix  [ ij  (a>)  ] 

*i 

- radius  vector  to  array  elements  i 

r 

*ij 

- gage  pair  position  vector 

S(w), 

(u)  - point  spectrum,  cross  spectra 

vi 
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T - length  of  time  record 

Wi(0Q;uo)  - weighting  coefficients  associated  with  beam  pattern 

W(t)  - taper  function  in  time  domain 

W(u>)  - Fourier  transform  of  W(t)  (window  function) 

x,y  - coordinates  in  physical  space 

xc  - position  vector  of  ith  sensor 

Ap  - pin  spacing  width  on  pin  type  wave  detector 

rij^Ct)  - time  history  associated  with  ith  wave  sensor 

0 - azimuthal,  angle  coordinate  associated  with  k 

X - waveleng ^h 

0 - standard  deviation 

4>ij  - argument  of  sensor  pair  (i,J)  vector 

a)  - frequency  in  radians /second 


I.  INTRODUCTION 


Oceanographic  literature  of  the  past  twenty  years  provides  a wealth 
of  descriptions  and  techniques  for  the  mathematical  modeling  of  ocean 
waves.  Since  the  acceptance  of  stochastic  processes  as  a viable  tool 
for  the  behavioral  prediction  of  surface  waves,  considerable  effort  has 
been  expended  in  an  attempt  to  measure  and  predict  the  structure  of  the 
wave  spectrum  associated  with  given  wind  speeds.  The  spectral  description 
for  time  histories  describing  the  surface  elevation  at  an  arbitrary  point 
are  fairly  well  developed  and  are  widely  used,  although  the  parameterized 
model  in  vogue  changes  rapidly  (cf  Hasselman,  1973  and  1976;  Pierson,  1976). 
However,  the  issue  as  to  how  the  spectrum  is  distributed  over  direction  is 
by  no  means  well  understood.  Current  knowledge  of  the  angular  spreading 
of  wave  spectra  is  minimal  and  wave  forecasting  awaits  a sufficient  data 
base  upon  which  to  build  satisfactory  directional  models. 

Essential  to  the  establishment  of  such  a data  base  are,  obviously, 
spectral  estimation  techniques  of  sufficient  resolution.  With  the  advent 
of  the  Fast-Fourier-Tranaform  (FFT)  and  data-adaptlve  spectral  analysis 
techniques  (e.g.,  the  MLM) , vast  quantities  of  raw  recorded  data  may  be 
handled  efficiently  by  cot . iters  such  that  temporal  processes  are  readily 
digested  and  spectrally  analysed.  Directional  spectra  naturally  involve 
spatial  transforms  and  spatial  cross  correlations  from  observations,  typic- 
ally over  a small  portion  of  the  wave  field.  Relatively  few  observation 
or  recording  points  can  cause  poor  angular  resolution  of  the  spectral  es- 
timate. This  may  be  observed  in  the  following  examples. 

Directional  wave  spectra  estimates  were  determined  from  measurements 
obtained  in  San  Francisco  Bay  and  extensively  reported  in  Chou,  £t  al, 

0974).  Two  contour  plots  of  wave  energy  density  spectra  S(u), 0)  are  shown 
in  Figures  1-1  and  1-2  where  frequency  increases  radially  and  wave  dir- 
ection corresponds  to  the  angular  variable.  Though  somewhat  difficult  to 
see  in  these  plots,  both  examples  suggest  that  there  was  a considerable  a- 
mount  of  directional  spreading  of  the  wave  energy.  On  the  other  hand,  ob- 
servers at  the  time  were  impressed  by  the  relative  narrowness  of  the  wave 
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distribution  and  the  lack  of  spectral  energy  from  direction  other  than 

v. 

the  principle  wind  wave  sources.  Some  of  the  ship  motion  predictions 
using  the  measured  directional  spectra  were  improved  if  the  spectra  were 
truncated  in  direction  so  as  to  Include  only  the  most  energetic  parts  of 
the  spreading  function.  Figure  1-2  is  a particularly  dramatic  example  of 
"cross  seas,"  two  wind  wave  systems  radiating  into  the  measurement  area 
at  right  angles.  The  contours  connecting  the  two  principle  peaks  are  most 
likely  spurious  and  are  due  to  the  analysis  of  the  particularly  small  ar- 
ray of  wave  sensors  used  inthe  experiment*.  Improved  resolution  can  be 
obtained  by  increasing  the  array  size,  though  often  at  considerable  cost. 
Since  the  spectral  analysis  was  done  using  the  data  adaptive  Maximum  Like- 
lihood Method  (MLM) , an  exceptionally  high  resolution  technique,  it  is 
possible  that  the  field  observations  and  ship  motions  theory  are  equally 
to  blame  for  the  apparent  lack  of  experimental  and  theoretical  correlation. 
An  analysis  of  the  resolution  potential  of  small  arrays  and  the  (nonlinear) 
MLM  spectral  estimation  techniques  appears  to  be  in  order  and  is  the  sub- 
ject of  this  report. 

Chapter  II  describes  some  of  the  wave  measuring  systems  and  array 
geometries  in  use.  Chapter  III  reviews  the  normal  recording  and  spectral 
analysis  steps  in  array  data  preparation  and  Chapter  IV  discusses  array 
resolution  in  terms  of  the  beam  pattern.  Chapters  V,  VI,  and  VII  provide 
the  wave  spectral  simulation  and  the  array  sensitivity  analysis. 
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II.  MEASURING  SYSTEMS 


In  order  to  determine  the  wave  number  spectrum  of  the  sea  surface,  it 
is  necessary  to  obtain  a measure  of  the  wave  amplitude,  or  related  varia- 
ble, over  a region  of  space.  Photographic  and  optical  reflection  tech- 
niques are  attractive  in  principle,  but  clumsy  to  work  with  in  practice 
(cf.  Cote,  et  al , 1960  and  Pierson  ,1976) . Many  of  these  techniques  are 

only  appropriate  for  the  short-gravity-wave  and  capillary-wave  regions 
(Cuthart,  et  al,  1975).  Land  and  satellite  based  wave  measuring  systems 
at  present  all  suffer  from  insufficient  resolution.  Only  fixed  spatial 
arrays  of  wave  sensors  are  repeatedly  used  with  success  in  the  spatial 
measurement  of  wind  generated  gravity  waves.  For  this  reason,  scanning 
techniques  and  the  combination  of  fixed  and  scanning  sensors  (cf.  Mollo- 
Christensen  & Koziol,  1972)  will  not  be  considered  in  this  report.  They 
not  only  require  special  treatment  in  the  interpretation  of  their  spectral 
analysis  due  to  Doppler  shift  effects,  but  they  also  require  special  con- 
sideration in  the  construction  of  the  covariance  matrix  described  in 
Chapter  IV. 

The  communications  and  slesmic  literature  is  replete  with  analyses 
and  optimum  design  procedures  for  large  arrays,  i.e.,  arrays  that  could 
contain  from  tens  to  hundreds  of  elements  (Capon,  1969).  However,  water 
wave  arrays  have  not  exceeded  six  elements,  with  three  and  four  sensors 
being  the  most  common.  These  arrays  are  extremely  small  and  certainly  do 
not  permit  approximate  design  procedures  based  on  a continuous  distri- 
bution of  elements.  A typical  array  configuration  is  shown  in  Figure  II-l 
along  with  the  wave  heading  angle  definition. 

The  most  common  wave  sensors  are  wave  staffs,  bottom  mounted  pressure 
transducers,  and  tethered  buoys.  The  wave  staffs  use  the  sea  water  as  an 
electrical  short  circuit  and  are  either  of  the  resistive  or  capacitive 
type.  The  former  allows  the  sea  water  to  modulate  the  resistance  of  a 
bare  wire  or  short-out  successive  pins  (step  type).  The  latter  is  de- 
signed to  measure  the  capacitance  between  a dielectrically  clad  wire  and 
the  water.  Pressure  transducers  are  attractive  because,  unlike  the  wave 


Fig.  II-l  Definition  Sketch:  Array  Pattern  with  four  elements 

staffs,  they  do  not  require  a special  (possibly  reflecting)  structure  on 
which  to  mount  them  and  are  relatively  safe  from  the  standpoint  of  vessel 
operations.  On  the  other  hand,  they  are  more  difficult  to  locate  pre- 
cisely and  have  poor  high  frequency  response.  Tethered  buoys  have  advan- 
tages in  deep  water  since  they  can  be  placed  where  a fixed  structure  would 
not  be  practical.  However,  extracting  the  wave  time  history  from  the 
buoy  motions  and  establishing  the  relative  position  of  the  members  of  the 
buoy  array  are  likely  to  be  difficult  tasks  at  best.  There  are  a large 
number  of  review  papers  describing  wave  measuring  and  monitoring  systems. 
Relatively  current  bibliographies  can  be  found  In  Pierson  (1976),  Guthart, 
et  al  (1976),  and  Panlcker  (1974).  A list  of  regularly  used  wave  measuring 
devices  is  given  by  DTaper  and  Fortnum  (1974). 

Numerous  array  geometries  have  been  used  and  it  Is  appropriate  to 
cite  a few  of  the  larger  and  more  recent  examples.  Bennett  (1964)  employed 
a six-gage  array  of  pressure  gages  in  the  Gulf  of  Mexico.  A five-gage 
array  was  installed  by  the  Corps  of  Engineers  off  the  California  coast 
(Borgman  and  Panicker,  1970)  and  the  configuration  is  shown  in  Figure  A-l 
(f)  of  Appendix  A.  More  recently,  a four-gage  array  of  wave  staffs  was 
used  in  San  Francisco  Bay  (Haddara,  et.  al,  1972)  and  the  resulting  direc- 
tional spectra  reported  by  Chou,  et  al  (1974).  Deep  Oil  Technology  re- 
cently completed  their  large  scale  Sea  Test  Program  of  a Tension  Leg  Plat- 
form off  of  Catalina  Island,  California.  A six-gage  array  of  wave  staffs, 
similar  to  the  one  sketched  In  Figure  A-l(e) , was  used  and  the  data  re- 
ported to  NOAA's  Environmental  Data  Service  (NODC,  1976).  Except  for  the 
last  two  test  programs,  the  analysis  of  the  data  obtained  from  the  arrays 
reported  to  date  have  been  analyzed  only  for  the  principle  wave  direction 
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or  with  "non-adaptlve"  spectral  techniques.  Capon  (1969)  has  shown  that 
potentially  far  greater  resolution  Is  possible  with  the  data-adaptlve 
maximum  likelihood  method.  This  report  will  explore  the  potential  im- 
provement in  resolution  when  the  MLM  Is  applied  to  data  obtained  from 
the  small  arrays  typical  in  ocean  wave  measurements. 
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III.  SPECTRAL  ANALYSIS 


1 . Temporal  Spectral  Estimation 

Computational  schemes  for  estimating  the  power  spectral  density  of 
a stationary,  random,  time-varying  recorded  signal  have  been  developed  in 
the  last  quarter  century.  Initially,  the  basic  procedure  was  to  form  an 
estimate  of  the  temporal  correlation  function  which  was  directly  transformed 
to  a "frequency-domain"  to  provide  an  estimate  of  the  spectral  density. 
This,  more  or  less,  is  the  methodology  advanced  by  Blackman  and  Tukey 
(1959).  In  1965,  Cooley  and  Tukey  Introduced  the  computationally  superior 
Fast  Fourier  Transform  (FFT)  technique  to  obtain  spectral  estimates.  This 
method  employs  the  product  of  the  transformed  time  series  rather  than  the 
conventional  transform  of  a synthesized  correlation  function.  Both  of  the 
above  methods  will  be  generally  referred  to  as  direct  spectral  estimators. 
For  example,  suppose  an  eetlmate  of  the  cross-spectrum  of  two  measured 

A A 

time  histories n ]^(t) , n 1*  desired  and  these  records  are  each  of  length 
T.  Let  A^ujjT),  the 

A 

Here,  the  ^(t)  and  A 
taper  function.  The 

This  estimate  has  the 


direct  Fourier  transform,  be  defined  as: 

'Qt***^-  a=  (3.i) 

^(uit)  are  random  variables  and  W(t)  is  a suitable 
(FFT)  spectral  estimate  is  then  given  by: 

— — — mm—  ■■  mem  s 

j <3*2) 

property  that: 
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where  is  the  true  cross-spectrum. 

Tapering  is  inherent  to  all  spectral  estimation  procedures  and  is 
applied  for  numerical  stability  of  the  spectral  estimate.  The  taper 
function  is  applied  directly  to  the  estimated  autocorrelation  function 
(Blackman-Tukey  method)  or  the  time  series  (FFT)  method  before  any  Fourier 
transforms  are  calculated.  The  shape  and  duration  of  the  taper  function 
are  determined  a priori  through  consideration  and  trade-offs  of  desired 
spectral  resolution  and  statistical  variance.  The  window  function  Is  the 
Fourier  transform  of  the  taper  function  and  relates  the  average  estimated 
spectrum  to  the  true  spectrum.  Proper  taper/wlndow  function  design  Is 
difficult  in  that  high  frequency  resolution  tends  to  be  achieved  at  the 
expense  of  statistical  reliability,  and  for  limited  amounts  of  data,  this 
phenomenon  is  particularly  troublesome  and  statistical  reliability  is 
particularly  poor.  The  dominant  characteristic  of  direct  or  conventional 
algorithms  is  that  the  taper  or  window  function  is  prescribed  and  does  not 
depend  on  the  true  spectrum  being  estimated, cf.  Otnes  and  Enochson  (1972). 

Within  the  last  ten  years  powerful  data  adaptive  methods  have  developed 
to  alleviate  the  windowing  problem.  These  methods  are  generally  termed 
"nonlinear"  vs.  the  direct  or  linear  methods  previously  mentioned.  The 
window  associated  with  these  methods  depends  on  the  true  spectrum  as  well 
as  the  power  at  the  particular  frequency  being  estimated.  It  adapts  in  an 
optimal  way  so  as  to  be  least  "disturbed"  by  noise  or  signal  power  in 
adjacent  frequency  bands.  Two  different  nonlinear  methods  are  currently 
available:  1)  the  maximum  likelihood  method  (MLM)  is  designed  to  minimize 

the  error,  in  a least  squares  sense,  of  the  spectral  estimate,  and  2)  the 
maximum  entropy  method  (MEM)  which  though  not  yet  state-of-the-art  has  the 
highest  potential  resolution  capabilities  (Lacoss,  1971).  The  MEM  is  some- 
what more  difficult  to  compute  and  the  spectrum  is  not  interpretable  in  the 
same  way  as  the  usual  power  density  spectrum.  It  is  therefore  not  as 
useful  a tool  for  engineering  analysis. 

The  MEM  method  is  still  being  refined  and  is  not  in  general  use.  Con- 
sequently, this  study  will  emphasize  the  direct  and  MLM  techniques. 


2 . Spatial  Spectral  Estimation 

The  meafurement  of  the  frequency-wave  number  spectrum  In  a multidi- 
rectional seaway  by  means  of  static  arrays  Is  done  by  employing  the  spatial 
analogs  of  the  temporal  estimation  schemes  discussed  previously.  The 
window  functions  are  now  addressed  to  wave  number  space  (magnitude  and 
direction).  The  directional  window  function  is  called  the  array  beam  pat- 
tern. The  beam  pattern  associated  with  the  direct  or  conventional  fre- 
quency-wave number  estimation  is  called  the  "natural"  beam  pattern.  As  in 
the  pure  time  series  analog,  the  natural  beam  pattern  does  not  adapt  to  the 
wavenumber  spectrum  being  estimated,  but  depends  entirely  on  the  geometric 
configuration  of  the  array.  Similarly,  an  "effective"  beam  pattern  is 
associated  with  the  adaptive  or  nonlinear  methods.  This  directional  win- 
dow adjusts  to  minimize  Interference  from  directional  components  in  neigh- 
borhoods adjacent  to  the  direction  of  immediate  interest. 

The  formulation  of  the  direct  and  MLM  spectral  estimates  are  given  as 
follows.  Covariance  matrix  elements  arc  defined  as: 


Cj.to'*  - Y (3.3) 


where  A1(u;T)  is  the  Fourier  transform  of  the  measured  wave  amplitude  n^(t) 
as  in  Equation  (5.1).  The  matrix  elements  are  generally  normalized  using 
the  deviations  oi(u>)  -t/sn"  (u)  in  order  to  reduce  the  effect  of  calibration 
errors,  i.e., 


Zn  m 

(SJWSG'.Cw') 


The  direct  frequency-wavenumber  spectral  estimate  Is  given  by 

s tojto  = e.  wiktoVr'i 

4«<  \*\  l\ 


(3. A) 


whereas  the  maximum  likelihood  estimate  of  the  spectrum  is  (Capon,  1969): 
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(3.5) 


i 

where (oj)  are  elements  of  the  inverse  of  the  matrix  II  Ci-CosA . 

J 1 /V 

For  our  purposes  it  is  assumed  that  the  covariance  elements  C1j 


are  reliable  estimates  since,  in  practice,  T,  the  record  length,  is 
typically  large. 


3.  Recording  and  Analysis  Procedure 

The  measurement  and  analysis  procedure  typically  employed  for  wave 
studies  is  shown  in  Figure  III-l.  Such  a format  is  implicit  in  the 
simulation  and  analysis  steps  discussed  below.  Detailed  considerations  in 
the  measurement  and  temporal  spectrum  analysis  steps  are  given  by  Bendat 
and  Piersol  (1971).  The  problem  of  spectral  estimate  variance  reduction, 
often  called  smoothing,  is  discussed  by  Jenkins  and  Watts  (1969)  and  Otnes 
and  Enochson  (1972). 

It  is  assumed  that  the  requirements  for  good,  normalized  estimates  of 
the  covariance  matrix  can  be  met.  The  remainder  of  this  report 

will  address  the  problem  of  directional  spectrum  resolution,  particularly 
for  the  MLM  applied  to  array  analysis. 
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FIGURE  III-1  Recording  ard  Analysis  Procedure 
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IV.  MATHEMATICAL  DESCRIPTION  OF  THE  BEAM  PATTERN 


The  array  is  modeled  as  a static  aggregation  of  the  elements  attached 
to  a fixed  rigid  reference  frame.  Each  element  in  turn  is  imagined  to  be 
a point  sensor,  e.g.,  in  the  simplest  case  a recording  device  which  re- 
sponds linearly  to  wave  elevation. 


I 


Wave  Elevation 


in 


ith 

SENSOR 

ELEMENT 


Voltage 


out 


FIGURE  IV-1 


The  recorded  voltages  from  each  sensor  element  are  subsequently  col- 
lected and  jointly  processed  by  digital  computer.  The  resolution  capabil- 
ity of  any  spatial  processing  scheme  may  be  inferred  by  its  characteristic 
beam  pattern  or  "directional  window."  The  quantity  used  here  to  denote  the 
beam  pattern  of  an  array  is 

BC0;  co  | 0,;coo] 


where  0 is  the  azimuth  or  spatial  variable ; U)  Q is  thecircular  frequency  at 
which  the  spectrum  is  being  estimated;  and  0q  is  the  direction  at  which 
the  spectrum  is  being  estimated. 

The  estimated  value  of  the  spectrum  is  related  to  the  true  value  of 
the  spectrum  explicitly  by: 

$<e.;co.^  =li  S(S,^BL9;<al4|«.ldedw 


mmm 
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where  S(0,m)  is  the  true  frequency-wavenumber  spectrum,  S(0o,o^)  is  its 
estimate  at  and  | denotes  the  expected  value  operator. 

Roughly  speaking,  the  most  desirable  beam  pattern  is  that  which  is  "smooth" 
about  the  estimation  point  (®0^0)  and  negligible  elsewhere,  i.e.,  in  the 
limit  a vanishingly  small  slit  centered  about  the  estimation  point.  As 
indicated  by  Expression  (4.1),  the  true  spectrum  is  "leaking"  through  the 
spatial  window  offered  by  the  beam  pattern  in  yielding  an  estimate  of 
S(eou,o).  The  exact  structure  of  the  beam  pattern  depends  entirely  on  the 
spatial  configuration  of  the  array  and  the  nature  of  the  processing  involved. 

In  general,  a (discrete  array)  beam  pattern  may  be  written  as. 


iM  M 1 1 


(4.2) 


where  N «*  number  of  sensor  elements 

r^  ■ Xj  - x^  • gage  pair  position  vector 
k(<u)  • the  wavenumber  vector 


- weighting  coefficients,  in  general  modified  for  each 


estimation  point  (0  , u)  ) • 
o o 


It  is  customary  to  define  the  weights  W^(0q,  Wq)  30  that: 

* 1. 


(4.3) 


The  simplest  form  of  (4.2)  is  the  so-called  "natural"  beam  pattern.  In 
this  case. 


kue.iw.3  - jj-  elfe“*V2; 


(4.4) 


where  ^(w)  “ ko[cos0o,sin0oJ  and  kn  is  obtained  from  the  dispersion  re- 


lation •K . The  natural  beam  pattern  becomes: 
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-izie1'  - ‘l  (4.5) 

N iM  j*l 

Thus,  the  natural  beam  pattern's  characteristics  are  determined  entirely 
by  the  location  and  number  of  its  sensors,  Mji'V"  A typical  example 
is  an  array  of  two  concentric  triangles  shown  in  Figure  IV-2. 


WAVE 

HEIGHT 
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Figure  IV-2 

The  array  is  more  or  less  symmetric  with  rings  of  three  elements  each  at 
radii  r^,  r2  respectively.  The  natural  beam  pattern  with  id  held  fixed 
is  given  approximately  by: 

s jL  + J#(£refe  3in(^*)'))2'  (4.6) 

where  Jq(  ) is  the  0th  order  Bessel  functions  and  -tank  fek. 

This  expression,  containing  the  lead  terms  describing  the  structure  of  the 
array,  exhibits  the  typical  characteristics  of  sparse  arrays.  In  parti- 
cular, as  rj , r^-¥-oo  , the  main  beam  narrows  in  the  neighborhood  of  0“6Q. 
which  is  desirable.  In  fact,  it  may  be  shown  that  the  "half  power"  beam- 
width  of  the  major  lobe  of  the  natural  pattern  of  any  array  is  given  by: 
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(4.7) 


Note  that  Eq.  (4.7)  Implies  that  the  beamwidth  is  inversely  proportional 
to  the  root-mean-square  distance  of  the  sensor  elements  from  the  center  of 
the  array.  However,  as  the  RMS  distance  becomes  larger,  sidelobe  levels 
become  more  significant  which  in  turn  degrades  the  array  performanfce 
through  leakage.  (This  is  not  evident  in  Expression  (4.6),  however,  as 
the  secondary  lobe  structure  has  been  suppressed  for  compactness).  The 
natural  beam  pattern  is  inherent  to  the  straightforward  direct  spec  i ^1 
estimation  techniques  and  the  design  problem  is  the  trade-off  between  beam- 
width  and  sidelobe  level. 

A concise,  introductory  treatment  of  the  (natural)  resolution  capa- 
bilities of  static  arrays  was  presented  by  Barber  (1961)  some  years  ago 
(a  synopsis  of  which  was  included  in  the  work  of  Barnett,  1969  ).  Some 
simple  arrays  that  are  examined  in  order  to  assess  their  resolving  power 
are  the  line  array,  the  triangular  array,  and  the  star  array.  As  non- 
linear spectral  estimation  methods  were  not  known  at  the  time,  the  under- 
lying framework  of  the  discussion  appeals  to  direct  methods  of  spectral 
estimation.  The  two  figures  shown  in  Figure  IV-3  are  reproduced  from 
Barber  (1961).  They  represent  contour  plots  of  the  beam  patterns  for  a 
triangular  and  a four-element  array  respectively.  The  coordinate  space 
corresponds  to  wavenumber  space  (k  - space). 

The  arrays  are  "steered"  to  a particular  wavenumber  k0«[l0,m0] 

The  directional  aspects  of  the  array(s)  may  be  ascertained  from  the  inter- 
section of  the  locus  of  points  on  the  circle  of  radius  |k^k0  with  the 
contours  of  the  beam  pattern  as  suggested  in  Figure  IV-3(e).  A more  de- 
tailed discussion  is  given  by  Barber. 

For  the  more  sophisticated  adaptive  techniques,  the  natural  beam 
pattern  does  not  apply.  One  must  Instead  introduce  the  effective  beam 
pattern.  The  MLM  technique  in  particular  weights  the  information  from 
each  sensor  in  an  optimal  manner  depending  upon  the  statistics 
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<>f  tlu*  wave  field  Itself  (l.e.,  "adaptive").  The  effective  beam  pattern 
for  the  maximum  likelihood  method  may  be  written  In  the  form: 


B TJ2  l4*.(9;co\0ej6O^  e 

MLM  ^ 4 


(4.8) 


A comparison  of  Eqs.  (4.5)  and  (4.8)  reveals  the  additional  weighting 
coefficients  in  the  latter.  As  the  functional  notation  implies, 

these  coefficients  vary  for  each  choice  of  and  are  unique  for  a 

given  array  geometry,  sensor  noise  and  wave  field.  The  result  is  that  the 
MLM  method  (as  well  as  other  adaptive  schemes)  can  potentially  resolve  a 
wave  spectrum  with  greater  overall  precision  than  the  traditional  direct 
spectral  analysis  techniques  which  rely,  in  effect,  on  the  natural  beam 
pattern  of  the  array. 

The  effectiveness  of  these  data  adaptive  techniques  can  be  demon- 
strated by  simply  plotting  representative  beam  patterns  generated  for  a 
particular  wave  field.  Consider,  for  simplicity,  a directional  spectrum 
of  the  form: 


S(co^  = 5(60^)  FC^ 


where  F(0)  is  the  "normalized  directional  spectrum  spreading  function," 
and  the  three  element  array  pictured  in  Appendix  A,  Figure  A-l(a).  The 
signals  received  at  each  sensor  of  the  array  will  be  assumed  to  be  noise 
free,  i.e.,  there  is  no  additional  energy  or  ambiguity  associated  with  the 
time  series  to  be  processed.  This  is  never  true  in  practice,  cf.  Chapter 
IV,  but  the  resulting  beam  patterns  are  a particularly  dramatic  example  of 
the  power  of  adaptive  techniques. 

Let  us  now  consider  what  properties  we  would  wish  the  beam  pattern  of 
an  array  to  possess.  Since  the  true  spectrum  is  multiplied  by  the  beam, 
and  integrated  over  direction,  we  would  like  the  beam  to  be  unity  in  the 
direction  of  interest,  say  ©o,  and  suppressive  at  any  other  direction. 

This,  of  course,  would  be  accomplished  if  the  pattern  was  zero  everywhere 
but  the  direction  of  interest,  i.e.,  the  "look"  or  target  angle.  The 


-18- 


t 


'X. 

3 


j 


t -4 

CO 


1 ;it ter  is  ;i  very  severe  mathematical  constraint  and  is  unnecessary  for  this 
|i.i rt  lew I .ir  wave  field.  The  beam  need  only  be  peaked  in  the  direction  of 
interest  and  minimal  wherever  else  there  is  energy  in  the  input  wave  field. 

o ° 

li  no  energy  exists  in  the  range  -90<Q<  90  for  example,  the  beam  may  take 
on  arbitrary  values.  The  MLM  beam  pattern  can  be  shown  to  be  of  the  form: 


Z El (S')  Q <<4  E„(« 
B[0;tol6o)a>n  = = — 


V 

m.n 


where  E^(0)  is  is  the  radius  to  the  array  elements,  the 

prime  denotes  complex  conjugate,  and  0^(0))  is  the  inverse  of  the  covari- 
ance matrix  l\(Baggeror,  1976). 

The  following  figures  show  selected  MLM  beam  patterns  for  Geometry 
i-  1,  i.e.,  the  three  element  array  shown  in  Appendix  A,  Figure  A-l(a).  The 
pattern  for  look  angle  0O  multiplies  the  input  spectrum  and  the  result  is 
integrated,  as  in  Eq.(4.1),  to  obtain  the  estimate  of  the  spectrum  in 
the  direction  0 = 0 . The  synthesis  techniques  employed  in  these  numer- 
ical examples  are  discussed  in  Chapter  V.  Figure  IV-4  shows  the  beam  pat- 
tern at  four  "look  angles"  = 0°,  135°,  150°,  and  180°)  for  the  case  of 
a plane  wave  of  frequency  w = 0.2  rad.  and  incident  angle  0 - 150°.  The 
beam  patterns  are  plotted  in  decibels,  i.e.,  10  log^QB[0,(J,  on  a linear 
scale  every  five  degrees  with  the  origin  as  -15db.  A small  arrow  just  out- 
side the  unit  circle  (Odb)  points  (approximately)  to  the  look  angle  0 . 

o 

Also  plotted  are  the  spreading  functions  or  normalized  directional  spectra 
for  the  input,  the  direct  analysis,  [Eq.(3.Q]and  the  MLM  analysis 
f Eq . ( 3. 5 ) | . The  input  here  is  a spatial  delta  function  and  is  represented 
by  a triangle  of  unit  area  and  is  drawn  with  a solid  line.  The  direct 
spectral  estimate  is  shown  as  a series  of  small  triangles  and  is  nearly  a 
constant  for  this  example.  The  MLM  spectral  estimate  is  shown  as  a solid 
line  with  "plus"  signs. 

Note  that  the  beam  pattern  is  unity  (Odb)  in  the  0 direction  as 

o 

required  by  Eq.  (4.3).  Outside  of  0^,  the  beam  pattern  is  effectively 
null  in  the  Input  energy  direction  but  can  vary  wildly  otherwise.  These 
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large  variations  "allow"  the  main  lobe  to  ba  narrow  but  are  of  no  other 

importance.  Relatively  little  energy  la  "laaked"  into  the  aatimate  of 

S(0  ) from  directions  other  than  0 . For  0 ■ 150°,  the  beam  pattern  1« 

o o o 

approximately  unity  for  all  angles  0.  This  again  is  of  no  importance 
since  the  input  spectrum  Is  zero  outside  of  0 - 0^  and  their  product  does 
not  impact  on  the  Integral  of  Eq.  (4.1). 

This  is  a particularly  dramatic  example  of  the  power  of  the  MLM 
over  the  direct  method  of  array  processing.  Here,  the  MLM  effectively 
reproduces  the  input  spreading  function  while  the  direct  method  cannot  re- 
solve the  input  at  all.  Figure  IV-5  shows  the  beam  patterns  and  normalized 
spectra  for  a somewhat  shorter  wavelength  with  frequency  oj  * 0.6  radians. 
The  MLM  estimate  is  all  but  indistinguishable  from  the  input  while  the 
direct  approach  can,  at  best,  resolve  the  input  direction  without  any 
measure  of  its  narrowness. 

This  illustration  suggests  that  a three  element  array  might  perfectly 
resolve  all  but  the  longest  waves.  Unfortunately  this  is  not  the  case. 

A three  element  array  is  too  sparse,  i.e.  there  are  insufficient  terms  in 
the  Fourier  series  of  Eq.  (3.5),  to  describe  a spreading  function  with  a 
much  more  complicated  structure.  Figure  IV-6  shows  the  beam  patterns 
and  spectral  estimates  for  such  a case.  Here  F(0)  has  been  normalized  to 
unit  area  and  is  given  by 

[§  cosz&  4-  0.1  ( fcCUcT)  -V-  S(  15o°yfl  /n 

9o°<  9 < ZTo* 

o \9\<9o° 

The  plane  waves  at  100®  and  150®  each  contain  ten  percent  of  the  energy 
2 

under  the  cos  0 spreading  term.  The  delta  functions  1(0)  are  presented 

graphically  as  triangles  with  the  correct  area  at  the  plotted  5°  spacing 

of  the  spectral  input.  At  this  wavelength  the  direct  method  does  a very 
, 2 

poor  Job  of  resolving  the  cos  0,  let  alone  the  plane  waves.  The  MLM's 
resolution  is  much  better,  but  this  small  array  does  not  have  the  struc- 
ture to  isolate  the  plane  waves.  More  to  the  point,  real  world  measure- 
ments are  not  ideal  in  that  large  amounts  of  noise  are  typically  present. 


F(fl  = i 


These  extraneous  signals  greatly  reduce  the  MLM's  ability  to  adapt  to 
the  data  and  result  in  significant  reduction  in  spectral  resolution. 

In  order  to  analyze  the  effects  of  such  things  as  wavelength,  noise,  and 
array  geometry  on  the  resolving  ability  of  the  MLM,  a number  of  pertinent 
topics  dealing  with  the  simulation  of  real  data  are  discussed  in  Chapter  V. 
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V . Till-:  CO  VARIANCE  MATRIX.  RESOLUTION.  AND  NOISE 

Due  to  the  nonlinear  and  adaptive  character  of  the  MLM,  the  amount 
of  directional  spectral  "resolution"  one  might  expect  to  achieve  with  a 
given  system  is  difficult  to  quantify.  Resolution  depends  on  the  array 
size  and  geometry,  the  frequency  range,  and  the  amount  of  "noise"  intro- 
duced during  the  measurements  and  analysis.  Analytical  estimates  of  the 
importance  of  each  of  these  factors  are  in  general  not  available  since  the 
properties  of  the  inverse  cross-covariance  matrix  are  not  well  known. 

In  order  to  provide  guidance  for  the  design  of  array  systems,  it  will  be 
necessary  to  resort  to  numerical  simulations  of  sample  cross-covariance 
matrices  for  use  with  the  maximum  likelihood  analysis.  This  will  require 
the  selection  of  representative  directional  spectral  models,  corrected  for 
the  effects  of  noise,  and  the  taking  of  their  spatial  transforms. 

1 . Synthesis  of  the  Covariance  Function 

The  simplest  ocean  wave  model  is  an  infinitely  long  crested  wave  of 
length  , amplitude  A,  and  circular  frequency  oj  moving  in  the  © direction, 
e.g.,  n(t,x)  « A cos  (k  *x  -wt)  where  k * [kcos©,ksin©]  and  k(  (j)  obeys 
the  dispersion  relationship.  The  directional  spectrum  is  singular  in 
frequency  and  direction  and  a particularly  useful  model  for  testing  the 
limits  of  the  type  of  resolution.  A more  realistic  model  of  the  sea  sur- 
face is  the  multi-directional  random  wave  model.  The  directional  spectra 
referred  to  in  this  report  are  characterized  by  their  separability  into  a- 
zimuth  (©)  and  frequency  (co).  Although  this  choice  is  mathematically  re- 
strictive and  does  not  conform  to  the  spreading  models  suggested  by  Pier- 
son (]d76),  it  includes  a sufficiently  large  number  of  instructive  cases. 

By  hypothesis,  it  is  assumed  that  the  directional  spectrum  S(i),©)  may  be 
decomposed  as 


• SCcoh 

(5.1) 


where  the  "spreading  function"  F(6)  embodies  all  the  directionality  of 
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S(<^0)  in  the  azimuth  dependence  0 and  S(ui)  contains  all  the  frequency/ 
wavelength  character.  Hence  the  spreading  function  F(0),  by  definition, 
lias  the  property  that 

Zir 

J Rtf)  = i.  (5. 

o 

* 

The  spatial  covariance  matrix  C(o),r)  where  jr  is  the  spatial  lag,  is 
given  by  the  inverse  spatial  transform  of  the  spectrum 


C 


T? 7 


(5.3) 


Note  that  k . r » kr  cos(0  - <J>)  where  r = i*m-  arg(r) , and  k is  the 
wavenumber  magnitude  demanded  by  the  dispersion  relation.  Since  the 
interest  here  is  in  discrete  arrays,  the  covariance  function  need  be 
evaluated  only  at  the  array  element  pairs  r^  ~ — J * *,e** 


C 
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(5.4) 


and  is  referred  to  as  the  covariance  or  cross-covariance  matrix.  The 
temporal  spectrum  S(w)  is  arbitrary  and  may  be  taken  to  be  any  desired 
point  spectrum.  Since  the  concern  here  is  limited  to  the  directional 
resolution  of  the  array,  the  signal  power  can  be  normalized  and  the  de- 
sired spreading  function  selected  for  each  frequency.  Formally,  the 
spreading  function  F(0)  can  be  expanded  in  a Fourier  series  in  0,  i.e.. 


Rtf-  i o«< z‘"6 

n»-oo 


(5.5) 


where  1 an)  are  the  Fourier  coefficients  given  by 
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Using  the  identity 


(5.6) 


±£YiarcosC0'*H"%  = ain(^^jocfe0 


(5.7) 


the  covariance  function  becomes 


LA /*"!  ^ 


n--« 


(5.8) 


where  the  coefficients  {a^  are  given  by  Eq.(5.6).  This  form  for  C^ai.r), 
while  instructive,  is  not  necessarily  computationally  efficient  and  more 
direct  schemes  for  its  evaluation  will  be  employed  where  expedient. 

Three  types  of  spectral  spreading  functions  and  the  resulting  spatial  co- 
variance  matrices  are  discussed  in  Appendix  B,  namely  an  omnidirectional, 
a raised  cosine,  and  a delta  function  spectrum.  These  are  sketched  in 
Figure  V-l  and  are  intended  as  basic  examples  of  directional  spectra  for 
use  in  evaluating  the  resolution  capabilities  of  various  array  geometries. 
These  spreading  functions  are  by  no  means  exhaustive  though  various  com- 
binations of  them  yield  an  ample  variety  for  present  purposes. 


2 . Wavelength  Resolution  and  Matrix  Stability 

The  essential  requirement  for  use  of  the  maximum  likelihood  method 
of  spectral  estimation  is  the  existence  of  the  inverse  (i.e.,  nonsingular- 
ity) of  the  covariance  matrix  (^(w).  In  theory,  any  realizable  propa- 
gating wave  field  will  have  a non-singular  covariance  matrix.  Unfor- 
tunately, this  is  far  from  the  case  in  terms  of  numerical  stability.  In 
fact,  both  the  covariance  matrix  and  its  Inverse  rapidly  become  computa- 
tionally singular  as  the  array  size  to  wavelength  ratio  tends  to  zero. 

This  can  be  easily  seen  for  the  limiting  case  of  a unidirectional  wave 
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FIGURE  V-l  Spreading  Functions  for  Spatial  Covariance 
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train.  The  covariance  matrix  is  obtained  from  Eq . (5.3)  with  elements 

^CoS)  ; ler» 

c.  CoS)  = — <2*  ' 1 

L\  TL'rr  (5.9) 


where  k = 2r/A.  In  the  limit  as  X+m>  , the  C^j(o))  elements  all  become  unity 
and  the  matrix  and  its  inverse  are  therefore  singular.  In  practice, 
numerical  problems  arise  for  a noise-free  covariance  matrix  when  the 
wavelength  is  one  to  two  orders  of  magnitude  larger  than  the  typical  array 
dimension  r„.  Real  systems  are  normally  contaminated  with  a certain  a- 
mount  of  noise  which  usually  adds  a small  number  to  the  (real)  diagonal 

of  C.  (w) , as  will  be  shown  below.  Fortunately,  this  effectively  prevents 

1 5 -5 

the  matrix  from  becoming  singular  and  is  why  the  value  F(9)  x 10  is  used 

in  the  numerical  computations  when  the  noise  is  said  to  be  zero.  However, 
non-trivial  noise  levels  can  have  an  adverse  effect  on  the  directional  re- 
solution. 

At  the  other  end  of  the  spectrum,  i.e.,  short  wavelengths,  the  co- 
variance  matrix  is  more  highly  diagonalized  and  therefore  more  numerically 
stable.  Unfortunately,  finite  arrays  cannot  resolve  indefinitely  short 
wavelengths  due  to  spatial  aliasing.  Figures  V-2  through  V-5  again  show 
the  beam  patterns  and  directional  spreading  estimates  for  Geometry  #1,  the 
three  probe  array,  for  a series  of  wavelengths.  As  the  wavelength  becomes 
shorter,  the  MLM  beam  pattern  exhibits  more  structure  and  is  thus  able  to 
follow  more  closely  the  actual  spreading  function.  This  is  also  true  for 
the  direct  method.  Equation  (3.4),  but  there  is  an  increasing  number  of 
ambiguous  peaking.  In  the  final  figure,  the  frequency  is  2.8  rad.  and  the 
wavelength  is  only  about  half  the  average  radius  of  the  array  elements  and 
less  than  one  third  of  the  gage  pair  distances.  The  direct  spectral  es- 
timate has  completely  degenerated  and  the  MLM  has  shown  signs  of  numerical 
instabi 1 I ty . 

A consideration  of  the  beam  pattern  will  help  to  explain  the  degen- 
erate performance  of  the  array  at  the  shorter  wavelengths.  The  expression 
tor  tie  natural  beam  pattern  will  be  used  in  the  following  discussion, 
though  the  conclusions  will  later  be  seen  to  apply  to  the  MLM  beam  patterns 
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a»  well.  The  natural  beam  pattern  for  an  N element  array  is  given  by 
Equation  (4.5).  An  equivalent  form  of  the  summand  is 


“»(  0-9.-  •fcj') 


(5.10) 


where  A is  the  wavelength  and  are  the  polar  coordinate  dis- 

tances between  the  i-j  gage  pair.  A Fourier  series  expansion  for  the  real 
part  of  this  expression  may  be  written  as 


o A 2il  S J (5.11) 

where  the  Fourier  coefficients,  J L(2Vf7;/^  are  the  nth  order  Bessel 

Zk  ) 

functions.  For  long  wavelengths,  where  1,  the  series  converges 

very  rapidly  since  the  coefficients  fall  off  as 

j / Z.’KfrN  at  LL*<kttLl-  ^ kcOj  i Z ... 


Jz%r  * 


(5.12) 


However,  as  approaches  unity,  the  Bessel  functions  behave  as: 


!X\  - cost  - W-  j)  >oA., 


(5.13) 


The  series  coefficients  are  now  comparable  in  size  and  the  pattern  becomes 
highly  structured  in  angle  9 , i.e.,  the  pattern  oscillates  wildly  for 
small  Intervals  about  9*90—  . Thus,  the  beam  pattern  begins  to  degen- 

erate or  become  "hashy"  if  more  than  a few  of  the  interprobe  spacings  r.. 
become  comparable  to  a wavelength . This  effect  is  generally  referred  to 
as  spatial  aliasing.  The  interprobe  spacing  for  Geometry  #1,  for  example, 
is  approximately  87  feet.  In  Figure  V-2,  * 0.45  and  some  rippling  of 

the  direct  spectral  estimate  is  apparent.  In  Figure  V-3, ^»/^=l,55and  the 
disastrous  effects  of  aliasing  on  the  direct  method  are  evident. 
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i.  Noise  Contamination 

It  is  inevitable  that  the  wave  elevation  signal  will  be  contaminated 
by  "noise"  during  the  process  of  measurement,  recording,  and/or  processing. 
Here,  "noise"  is  taken  to  refer  to  all  aspects  of  the  measurement  a- 
n.ilysis  not  associated  with  the  propagating  wave  field  of  interest.  These 
would  include  such  things  as  inaccuracies  of  the  sensing  elements,  quanti- 
zation, round-off,  and  electronic  noise.  Signals  that  propagate  across  the 
array  but  are  not  associated  with  the  wave  field  that  is  being  measured 
are  referred  to  as  "propagating  noise."  Non-propagating  noise,  such  as 
60cps  electronic  noise,  is  often  of  very  high  frequency  compared  to  the 
wave  field.  Normally,  such  contamination  is  of  little  concern.  The  noise 
can  often  be  filtered  out  or  effectively  isolated  during  the  processing. 

A simple  increase  in  the  sample  size  will  usually  improve  the  statistical 
reliability  without  requiring  a redesign  of  the  entire  system.  Unfortun- 
ately, the  MLM  of  array  processing  is  particularly  sensitive  to  the  quantity 
of  noise.  A noise-to-signal  ratio  of  even  a few  percent  can  greatly  re- 
duce the  angular  resolution  of  the  small  arrays  in  question.  In  order  to 
estimate  the  resolving  capability  of  actual  arravs,  a noise  model  must  be 
included  in  the  simulation  and  analysis  process,  (cf  McDonough,  1972). 

There  are  numerous  potential  sources  of  signal  contamination,  some 
of  which  are  discussed  below.  The  two  types  of  noise  actually  addressed 
in  this  report  are  "sensor  quantization  or  non-propagating  noise"  and 
"propagating  noise."  Both  processes  are  modeled  as  independent  additions 
to  the  ambient  signal  field,  i.e., 

(5u) 

win  re  V7  is  the  measured  signal  and  n is  the  noise  component.  This  as- 
sumed independence  yields  simple  spectral  and  covariance  addition  of  the 
noise  components,  i.e.. 
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Sc^-)  = A-  J (5.15n) 

C'c«,c^  = C1Cu)JC;i')  V Cn(oJJQ.').  (5.15b) 

The  essential  problem  now  Is  to  obtain  an  estimate  of  the  slgnal-to-noise 
ratio  and  covariance  function  of  actual  array  systems. 

High  Frequency  Noise 

Wave  frequencies  above  the  bandwidth  of  Interest  and  high  frequency 
electronic  noise  can  pose  serious  aliasing  problems  if  they  are  not  effec- 
tively filtered  out  of  the  signal.  This  is  normally  done  prior  to  A/D 
conversion;  otherwise  the  filtering  must  be  done  digitally  with  a much  finer 
time  step  than  is  required  for  the  spectral  analysis.  Since  such  filtering 
is  normally  employed,  frequency  aliasing  will  not  be  discussed  further. 

The  problem  of  spatial  or  wavenumber  aliasing  will  be  addressed  below. 

Prop agating  Noise 

Propagating  noise  refers  to  any  signal  that  is  coherent  across  the 
array  yet  is  not  part  of  the  wave  field  of  interest.  The  major  source  of 
propagating  noise  for  a fixed  array  is  likely  to  be  the  reflections  from 
local  structures.  These  reflected  waves  represent  spatially  inhomogeneous 
signals  that  are  processed  like  any  other  propagating  signal.  The  inter- 
pretation of  the  measured  directional  spectrum  as  a spatially  stationary 
process  may  not  be  possible  if  there  is  a great  deal  of  local  interference. 
The  spectral  representation  of  locally  generated  (inhomogeneous)  wave  fields 
will  likely  be  highly  inaccurate  due  to  the  sparseness  of  the  arrays  in 
quest  ion . 

The  only  propagat tng  noise  component  of  the  synthesized  signal  will 
he  an  omnidirectional  term.  There  appears  to  be  no  evidence  that  might 
suggest  a significant  amount  of  naturally  occurring  wave  energy  distributed 
uniformly  in  azimuth.  The  directional  spreading  could,  of  course,  be  very 
broad  and  contain  numerous  humps;  nevertheless,  the  truly  omnidirectional 


J 


-14- 


component  is  assumed  to  be  small  and  essentially  uncorrelated  with  the 
principal  wave  field  of  interest.  As  shown  above,  a uniform  spreading 
function  F (9)  * 1/27T yields  a cross-covariance  term  of  the  form 


wiiere  N((,j)  is  the  noise  power  level  and  J (.)  is  the  zero  order  Bessel 

o 

function.  In  Chapter  VI  it  will  be  shown  that  omnidirectional  (propagating) 
noise  appears  to  be  treated  like  any  other  propagating  signal  by  the  MLM 

processing  scheme. 

Field  Noise 

The  desired  wave  field  signal  is  often  contaminated  with  random 
signals  due  to  spray,  wave  breaking,  sheltering  and  reflections  from  local 
structures,  water  runoff  from  the  sensors,  etc.  These  events  may  yield  a 
(noticeable)  contribution  over  a broad  band  of  frequencies  but  should  only 
be  weakly  correlated  across  the  array.  Their  effect  on  the  maximum  like- 
lihood processing  is  most  pronounced  when  the  wave  amplitudes  are  small  and 
there  is  a significant  noise-to-signal  ratio.  A priori  estimates  of  the 
magnitude  of  the  noise  power  are  difficult  to  make.  Most  wave  height  sen- 
sors peform  some  filtering  and  are  poorly  responsive  to  the  effects  of  spray 
and  wave  breaking.  Water  runoff  can  introduce  a lag  in  the  measured  signal 
and  distort  the  peaks.  One  would  not  expect  these  to  be  particularly  ser- 
ious problems,  however. 

On  the  other  hand,  contaminating  waves  generated  by  local  structures 
may  vary  significantly  in  amplitude  as  they  cross  the  array.  Local  re- 
flections and  sheltering  could  also  yield  signals  that  would  not  be  uni- 
formly coherent  across  the  array. 

f >i ta n t l /. i n£  Not se 

Quantizing  errors  can  arise  when  a continuous  variable  is  handled 
in  .1  discrete  form.  Round-off  errors  in  t he  course  of  digital  signal  pro- 
cessing on  a mini-computer  are  an  example.  In  most  situations  machine  round- 
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ot I errors  are  not  significant  or  can  be  reduced  by  going  to  higher  pre- 
cision. A more  relevant  example  is  the  quantizing  effects  of  a pin  type  wave 
detector  as  used  in  the  San  Francisco  Bay  experiments,  (Chou,  £t  al,  1974). 
This  type  of  wave  probe  is  able  to  resolve  wave  elevations  only  to  within 
the  pin  spacing, Ap,  on  the  staff.  This  anomaly  may  be  modeled  as  a type 
of  wideband  sensor  noise  which  effects  only  the  diagonal  of  the  covariance 
matrix  \\  For  portions  of  the  input  spectrum  with  component  am- 

plitudes several  times  greater  than  the  pin  spacing  width  "V  it  may  be 
^shown  that  for  a wave  sampling  rate  At,  an  effective  noise  power  Ap/lZAt 
is  added  to  the  main  diagonal,  C^Q^),  of  the  covariance  matrix.  This 
result  is  based  on  the  assumptions  that  1)  the  uncertainty  between  the  true 
and  measured  wave  elevation  is  uniformly  distributed  over  the  spacing  width 
Ap  2)  that  the  sampling  rate  is  on  the  order  of  the  "ttyquist"  rate  (i.e., 
the  reciprocal  of  twice  the  highest  significant  frequency  component  in  the 
spectrum),  and  3)  the  quantization  uncertainties  between  wave  probes  are 
uncorrelated.  Thus,  for  this  type  of  wave  probe  the  noise  matrix  is 

II  N \1  - dlcuji  ApV  12  A . 


An  expression  of  the  percentage  noise  is  given  by: 


ApV  12  At  , 

— E . \oo 

ft,  i3 


and  for  typical  values  of  say  Ap  = i"  , At  “(1/2  x 2.5  cps) and  very  small 
waves,  i.e.,  (ZXV*  $ S(<*»)ico  * oAlb  the  percentage  error  is  found  from  the 
above  expression  to  be  2.3%.  A well  designed  system  will  therefore  see  at 
most  a few  percent  quantizing  noise  for  the  smallest  waves  and  considerably 
less  for  the  design  situation. 
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VI.  SENSITIVITY  ANALYSIS 

■ 

Using  the  simulation  techniques  discussed  in  Chapter  V,  the  effects 
of  noise,  wavelength,  array  size,  etc.  can  now  be  discussed.  Numerous 
examples  have  been  studied  and  those  most  representative  or  illuminating 

are  presented. 

The  following  figures  include  plots  of  the  directional  spectral 
estimates  of  the  MLM  and  direct  methods  overlapping  the  input  spreading 
function.  In  addition,  they  depict  the  adaptive  beam  pattern  of  the  MLM 
at  four  distinct  look  angles.  Each  figure  is  associated  with  distinct 
geometry,  wavelength,  noise  figure,  and  input  spreading  function  (c f Chap- 
ter IV).  The  array  dimensions  and  frequencies  chosen  for  presentation  are 
appropriate  for  typical  ocean  wave  spectra,  but  can  be  easily  transferred 
to  other  scales  simply  by  normalizing  with  the  wavelength. 

1 . Noise,  Wavelength,  and  Aliasing 
a)  Instrument  Noise 

Nonpropagating  noise  is  a signal  that  is  uncorrelated  across  the 
array,  i.e.,  Cn(u),r^j)  “ in  Eq.  (5.15).  If  the  noise  originates  in  the 
sensors  and  is  independent  of  the  incident  wave  field  of  interest,  it  can 
be  modeled  as  simply  an  addition  to  the  diagonal  of  the  covariance  matrljc. 
The  following  figures  portray  the  capacity  of  a three  element  array  to  re- 
solve both  plane  waves  and  distributed  inputs  impinging  on  the  array  in  the 
presence  of  varying  instrument  noise  levels. 

Figure  VI-1  shows  the  dramatic  effect  of  one  percent  instrument  noise 
on  the  MLM  spectral  estimate  when  compared  with  Figure  IV-4.  This  corres- 
ponds to  an  exceptionally  long  wavelength;  fifty  times  the  radius  of  the 
array.  Figure  VI-2  is  Identical  except  that  the  wavelength  is  now  about 
equal  to  the  distance  between  the  array  elements.  The  direct  estimate  is 
completely  aliased  and  the  MLM  estimate  shows  some  aliasing  at  180°  from 
the  Input.  It  is  clear  that  the  effect  of  noise  is  to  smear  the  true 
spectrum  about  the  angle  of  incidence  and  increase  the  directional  ambi- 
guity for  the  short  wavelengths. 
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Figures  VI-3  through  VI-5  are  again  for  a single  plane  wave  input  at 
6 * 150°  and  various  wavelengths,  but  the  instrument  noise  level  is  now 
ten  percent  of  the  plane  wave  energy.  It  is  clear  that  some  care  must  be 
taken  in  estimating  the  input  spectral  width  on  the  basis  of  either  the  di- 
rect of  MLM  estimates  for  such  a small  array.  Figure  VI-4  shows  slight 
aliasing  of  the  MI  M estimate  and  Figure  VI-5  has  slightly  increased 
aliasing  over  the  one  percent  case  of  Figure  VI-2. 

Figure  VI-6  through  VI-10  are  examples  of  the  three  element  array’s 
resolution  for  a more  complicated  input.  Here,  the  spreading  function 
F(d)  Includes  two  plane  waves  and  a cosine  squared  distribution.  The 
plane  waves  are  each  ten  percent  of  the  energy  of  the  cosine  spreading 
function.  None  of  the  examples  exhibit  much  ability  to  resolve  the  plane 
waves  and  the  high  frequency  case  ( cu  * 1.5  rad)  has  degenerated  entirely. 
Hence,  the  high  frequency  limit  of  the  analysis  depends  not  only  on  the 
wavelength  and  array  geometry,  but  on  the  character  of  the  spreading  function. 

It  has  been  suggested  that  noise  levels  could  be  considerably  larger 
than  the  one  and  ten  percent  that  have  been  shown  so  far.  Consider,  for 
example,  a six  element  array  (Geometry  #3)  with  u>  - 0.8  radians  and  a cosine 
spreading  function  raised  to  the  flfieth  power,  i.e., 

FCO-^  * °C  \Q-t8o  \ £ <30. 


Figures  Vl-lOa,  10b,  and  10c  show  the  spectral  estimates  for  instrument 
noise  levels  of  50,  100,  and  200X  of  F(0).  In  each  case,  the  spectral 
peak  has  been  noticeably  lowered,  largely  due  to  an  apparent  increase  in  an 
omnidirectional  component.  If  the  omni-directional  component  is  removed, 
the  remaining  estimate  is  still  wider  than  the  input  spreading  function. 

For  the  200%  noise  level  in  particular,  a fit  to  the  modified  MLM  spectral 
estimate,  i.e.,  ~ requir®»  a cosine  raised  only  to  the 

eleventh  power,  rather  than  the  fifty  of  the  input.  The  beam  patterns  for 
these  examples  are  very  broad  and  are  the  cause  of  the  poor  resolution. 

Figure  1-2  is  the  MLM  estimate  of  "cross-sea"  data  discussed  in  Chap- 
ter I.  Two  wave  systems  were  generated  in  different  regions  and  radiated 
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into  the  measuring  area,  creating  a classic  "waffle  pattern."  Assuming 
the  resulting  system  at  f = .4  Hz  could  be  approximated  by  two  plane 
waves,  a noise  estimate  of  the  measured  data  can  be  obtained  by  simulation. 
The  array  configuration  corresponds  to  Geometry  II 2 when  scaled  by  a factor 
of  about  ten.  Figure  VI-11  shows  the  simulated  spectral  estimate  for  ten 
percent  instrument  noise.  Note  that  the  direct  method  is  unable  to  differ- 
entiate between  seas  separated  by  ninety  degrees!  Figure  VI-12  shows  the 
MLM  spectral  estimates  of  the  simulated  data  for  one  and  ten  percent  in- 
strument noise  with  the  MLM  spectrum  of  the  real  data  obtained  from  Figure 
1-2.  The  agreement  between  the  measured  spectrum  and  the  estimate  of  the 
noisy  simulation  permits  an  interpretation  of  the  noise  level  for  this  par- 
ticular case.  The  noise  content  could  be  expected  to  be  abnormally  high 
for  this  example,  since  the  wave  amplitudes  were  exceptionally  large, 
breaking,  and  confused.  Hopefully,  ten  percent  is  a practical  upper  bound 
for  the  quantity  of  noise  likely  to  be  encountered  under  normal  operating 
conditions . 
b)  Propagating  Noise 

The  only  propagating  noise  that  has  been  modeled  is  the  independent 
omni-directional  distribution  discussed  in  Chapter  V.  The  effect  of  an 
omni-directional  term  is  simply  to  raise  the  level  of  the  distribution  and 
is  handled  without  apparent  difficulty  by  the  MLM  analysis. 

2 . Geometry 

a)  Array  Size 

The  resolution  capabilities  of  an  array  increase  rapidly  as  the 
number  of  elements  N is  increased.  On  the  other  hand,  large  arrays  are 
expensive  to  construct  and  maintain  and  should  be  worth  the  added  cost  in 
resolution.  In  order  to  indicate  the  type  of  improvement  that  might  be 
expected  for  such  sparse  arrays,  a number  of  other  geometries  have  been 
investigated.  Figures  VI-3,  VI-13,  and  VI-14  display  the  spectral  estimates 
of  a three,  four,  and  six  element  array  (Geometries  #1,  112,  and  03  respec- 
tively) in  identical  wave  fields.  The  signals  are  all  subject  to  ten  per 
cent  instrument  noise  and  w * 0.5  radians.  The  improvement  in  the  MLM 
spectral  estimate  is  evident,  particularly  for  the  six  element  array. 
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This  is  primarily  due  to  the  increase  in  N rather  than  the  slightly 
larger  average  array  size  (cf.  Appendix  A). 

Figures  VI-7  and  VI-15  each  have  cosine  squared  distributions  with 
plane  wave  additions  with b) * * * * * * * * *  (l)  = 1.0  radians  and  ten  percent  instrument  noise. 
The  six  element  array  shows  a distinct  improvement  over  the  three  element 
array.  There  is  more  shape  in  the  neighborhood  of  the  peak  and  less  omni- 
directional noise  power.  The  main  lobe  of  the  beam  pattern  of  Geometry  II 3 
is  narrower  and  there  is  far  more  structure  than  for  Geometry  II 1. 

Vast  improvements  are  possible  if  the  array  size  is  increased  to  the 
ten  elements  of  Geometry  II 5.  This  configuration  is  a regular  one  with  the 
elements  spaced  evenly  on  two  rings.  The  array  geometry  statistics  are 
substantially  the  same  as  arrays  II 3 and  II 4 and  are  given  in  Appendix  A. 
Figures  VI-16  and  VI-17  are  for  the  same  distributed  input  considered 
above;  w = 1.0  radians,  and  the  one  and  ten  percent  instrument  noise 
levels.  The  low  noise  level  of  Figure  VI-16  finally  allows  the  plane  waves 
to  become  visible  whereas  they  were  only  hinted  at  in  the  six  element  case. 
The  larger  noise  level  of  Figure  VI-17  reduces  the  resolution  of  the  plane 
waves  but  is  not  seriously  contaminated  by  omni-directional  energy.  The 
effect  of  noise  is  clearly  evident  in  the  beam  patterns. 

b)  Element  Spacing 

The  effects  of  regular  vs.  irregular  array  spacing  have  been 

investigated  by  comparing  Geometry  II 3 with  Geometry  HU.  The  latter  is  a 

six  element  "double  ring"  array  whereas  the  former  is  a somewhat  irregular- 

ly spaced  version  that  was  actually  employed  in  the  full  scale  measurements 

reported  to  the  NODC  (1976).  Figure  VI-18  is  the  double  ring  array  with 

ten  percent  noise  and  a wavelength  about  twice  the  RMS  gage  pair  distance. 
There  is  a fair  amount  of  structure  evident  in  the  beam  pattern  and  in  the 

estimated  MLM  distribution.  The  corresponding  case  for  the  irregular 

geometry  is  shown  in  Figure  VI-15.  There  is  a slight  decrease  in  the  amount 
of  aliasing  at  the  expense  of  some  resolution.  These  differences  are 
evident  in  the  variations  of  the  beam  patterns.  In  general,  a regularly 

spaced  array  will  have  a beam  pattern  with  a narrow  main  lobe  but  rather 

pronounced  side  lobes.  If  the  array  spacing  is  perturbed  so  that  it  is  no 
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longer  periodic,  the  side  lobes  ure  typically  decreased  at  the  expense  of  a 
somewhat  broader  main  lobe.  This  effect  can  be  seen  by  comparing  the 
patterns  presented  In  Figures  VI-15  and  VI-18.  The  decrease  In  resolution 
of  Geometry  // 3 Is  therefore  due  to  the  wider  main  lobe.  Geometry  //4 , while 
having  a somewhat  superior  resolving  power  for  the  present  case,  does  ex- 
hibit a tendency  toward  aliasing  due  to  the  large  sidelobe  levels. 

It  is  not  surprising  that  the  spectral  estimates  are  quite  sensitive  / 

to  the  array  spacing  used  in  the  computations.  The  covariance  matrix  is 
constructed  without  regard  to  the  array  dimensions  and  it  is  only  in  the  j 

last  stage  of  the  analysis  tn»t  the  geometry  is  introduced.  An  error  in 
the  measured  dimensions  or  a mistake  in  the  numbering  system  that  in  ef- 
fect shuffles  one  or  more  array  elements  can  have  a disastrous  effect.  An 
example  has  been  constructed  by  spectral  analyzing  the  covariance  matrix 
generated  for  Geometry  // 4 with  the  array  of  Geometry  //3.  The  standard 
analysis  of  the  six  element,  double  ring  array  is  presented  in  Figure  VI-18. 

If  this  same  covariance  matrix  is  analyzed  with  the  irregular  spacing  of 
Geometry  //3  the  spectral  estimate  is  radically  altered  as  shown  by  Figure 
VI-19.  The  spectral  peak  is  slightly  shifted  and  is  considerably  narrower. 

If  the  investigator  was  expecting  a narrow  and  strongly  centered  distri- 
bution, such  an  error  in  the  computation  geometry  could  easily  go  undetected. 

3.  Summary 

Noise  contamination  of  the  signal  processing  clearly  has  adverse  effects 
on  the  resolution  capabilities  of  the  MLM  technique  for  any  array  config- 
uration. This  is  manifested  by  the  undesirable  broadening  of  the  main  lobe 
of  the  beam  pattern.  The  overall  effect  of  higher  noise  levels  on  the  MLM 
is  to  blur  and  disperse  the  fine  structure  of  the  directional  spectrum. 

The  fractional  decrease  in  performance  due  to  noise  is  highest  for  denser 
arrays. 

The  resolution  capabilities  of  an  array  increase  substantially  with 
the  number  of  elements  and  results  in  a narrowing  of  the  main  lobe.  This 
is  particularly  true  for  sparse  arrays.  The  spacing  between  elements  may 
be  regular  or  irregular  with  the  former  characterized  by  a narrower  main 

La  | 

lobe  but  somewhat  higher  side  lobes  than  the  latter. 
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Spectral  estimated  are  seriously  impaired  by  erroneous  numerical 
perturbations  of  ar'ay  topology.  A shuffling  or  permutation  in  the  order- 
ing of  the  array  elements  results  in  serious  distortion  of  the  spectral 
estimate.  This  may  be  an  aggravating  problem  in  the  analysis  of  array  data 
obtained  from  buoys,  especially  for  short  wavelengths. 

For  a given  noise  level  and  arrav  dimension  there  is  a limit  to  the 
wavelength,  roughly  equivalent  to  the  average  element  spacing,  below  which 
spurious  results  may  be  expected.  This  anomaly  is  analogous  to  the  prob- 
lem of  frequency  aliasing  inherent  to  sampling  rates  below  the  Nyquist  rate. 

Obviously,  many  variables  impact  on  array  performance  such  as  a speci- 
fic input  spectrum,  array  size,  etc.  One  is  generally  faced  with  competing 
designs  and  configurations  for  particular  applications.  A useful  design 
tool  to  aid  in  the  selection  of  a suitable  candidate  is  affored  by  a 
"figure  of  merit."  Such  a quantitative  index  of  performance  is  valuable  in 
assessing  the  combined  effect  of  wavelength,  noise  array  configuration  and 
directional  distribution  of  input  energy.  A convenient  figure  of  merit  will 
be  developed  in  Chapter  VII. 
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VII.  FIGURES  OF  MERIT 

The  overall  performance  of  the  array  is,  naturally,  a function  of  the 
intrinsic  geometry  and  sensor  characteristics  as  well  as  the  particulars  of 
the  ambient  wave  field.  Selecting  the  most  suitable  array  geometry  among  a 
variety  of  configurations  is  a difficult  task  indeed,  due  in  part  to  the 
nonlinear  nature  of  the  MLM.  While  no  optimization  technique  is  currently 
available  for  array  selection  (it  may  not  even  be  worth  the  effort!)  a 
simple  index  or  measure,  i.e.,  a figure  of  merit,  is  necessary  to  differ- 
entiate between  candidates. 

Initially,  the  number  of  elements  as  well  as  their  spacing  are  varia- 
bles which  must  be  chosen  in  a manner  consistent  with  the  expected  wave 
field.  Obviously,  the  wave  spectrum  is  not  known  a priori  (except  perhaps 
in  a gross  sense)  and  we  are  forced  to  choose  one  or  more  design  (dir- 
ectional) wave  elevation  spectra  and  assess,  via  analysis  and/or  simu- 
lation, the  array  performance.  Once  the  designer  has  chosen  his  design 
spectra,  he  is  faced  with  the  problem  of  rating  the  performance  of  the 
array  candidates.  This  is  most  easily  accomplished  by  simulation  where 
the  (simulated)  "estimated"  spectrum  is  usually  compared  with  the  input 
or  design  spectrum  over  all  directions  and  frequencies.  However,  it 
seems  apparent  that  a quantitative  performance  index  or  a "figure  of 
merit"  is  desirable  in  that  it  is  a more  efficient  scheme  than  visual 
inspection  of  reams  of  print-out.  The  figure  of  merit  must  be  a consis- 
tent estimate  of  performance  and  inevitably  must  reflect  the  average  re- 
solution of  the  array  for  a given  design  spectrum.  To  this  end,  consider 

the  frequency  variable  (to)  fixed  for  this  discussion.  Since  we  are 

o 

concerned  here  only  with  the  spatial  spreading,  let  B [0; n/ represent 
the  natural  beam  pattern  of  the  array  as  defined  in  Chapter  IV.  Let  F(0) 
represent  the  spreading  function  which  has  been  normalized  such  that 

I2.tr 
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In  order  to  fix  ideas,  consider  the  natural  beam  pattern.  In  this  case, 
Ul8;  (ii/Q  5 d)  1 is  ideally  relatively  narrow  or  peaked  about  the  look  angle 

C>  i) 

8 = 8 and  the  spectral  estimate  (eq.  7.1)  may  be  written: 

ir 

El §<«.,"$  * B[6;4 16 . (7.D 

o 

Furthermore,  for  this  natural  beam  pattern,  the  integral 

-ir 

may  be  interpreted  as  the  effective  beaniwidth  of  the  array  at  ft  - fte.  In 

VJ 

practice,  however,  MLM  beam  patterns  are  not  peaked  enough  and  hence  do 
not  always  satisfy  the  approximation  of  Eq . (7.1)  for  all  8^.  A figure 
of  merit  should  be  a measure  of  the  divergence  from  good  resolutions, 
hence  define  the  quantity  1(8^;^): 

IT 

l(0o;co^  * \ (7,3) 

-t r 


Some  consideration  of  Eq . (7.3)  leads  to  the  conclusion  that  the  narrower 

the  effective  beamwidth  about  8 « 8 the  smaller  the  value  of  1(8  ; m) . 

o C (D 

Now  the  fraction  of  energy  incident  in  a neighborhood  8^  of  a given  look 

angle  8=8  is  1(0  ; w).  It  follows  that  a "weighted"  measure  of  array 

o o o 

performance,  call  it  P(0j;u3j),  is: 

Pce0;cO  = Re.Vi(eo;^  (?.*) 

and  thus  the  performance  averaged  over  all  look  angles  046.4  2ir  is : 

Ztt  iy 


(7.5) 
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Wf  mny  define  /in  nrrny  gain  by  taking  the  ratio  (in  db's)  of  ?(*>  ) of 
tlie  given  array  to  that  of  a single  element  which  is  Inherently  omni- 
directional. For  a single  element  B[9;u)|6  ■ 1,  all  (0,  Gq)  , and  ((u,uio) 

from  Eqs . (7.3)  through  (7.5),  it  follows  that: 

[ P6so*V]  airUji4  2^  * (7 .6) 

elcmani 

Therefore,  the  array  gain  is  simply 

From  Eq.  (7.6)  it  follows,  for  example,  that  the  array  gain  of  a 
single  element  is  0 DB's  while  that  of  an  array  with  say  twice  the  average 
resolution  is  3 DBs.  The  array  gain  is  a measure  of  the  average  perfor- 
mance (i.e.,  averaged  over  all  look  angles  - Tf  £ 

We  choose,  therefore,  in  the  discussions  to  follow,  to  define  the 
figure  of  merit  (denoted  F.O.M.)  as  the  array  gain  in  decibels,  i.e. 

FTO.H  - 10  laj  (2/t rV(u>£).  a.D 

Though  the  ability  of  a given  array  to  resolve  directional  spectral 
features  at  various  wavelengths  (or  frequencies)  is  highly  geometry  depen- 
dent, resolution  trends  should  be  the  same  regardless  of  the  particular 
algorithm  (direct  or  adaptive)  implemented.  In  the  case  of  the  natural 
beam  pattern,  for  example,  expression  (4.8)  exhibited  the  inverse  re- 
lationship between  the  half-power  beamwidth  and  the  root-mean-square  dis- 
tance of  the  sensor  elements  from  the  arrays  geo-center. 

Figure  VII-1  shows  a plot  of  the  F.O.M.  vs.  normalized  wavelength  for 
four  arrays.  The  wavelength  has  been  scaled  by  the  average  interprobe 
distance  of  each  individual  array.  The  directional  spreading  of  the 
incident  wave  field  was  modeled  as  a mixture  of  two  delta  functions  and  a 
cosine-squared  function.  It  was  felt  that  this  directional  distribution 
was  a reasonable  combination  with  which  to  test  the  performance  of  the 
MI.M  algorithm  in  conjunction  with  the  various  array  geometries.  As 
'■xpected,  there  1b  a notable  improvement  in  performance  as  the  number  of 
sensors  Increases,  especially  in  that  region  where  the  wavelength  is  on 
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the  order  of  the  average  interprobe  distance.  In  particular,  when  the 
wavelength  is  about  equal  to  the  interprobe  distance  the  rate  of  increase 
of  the  F.O.M.  is  roughly  aaymptotic  to  .45  db's  per  sensor.  As  the 
wavelength  grows  large  compared  with  the  interprobe  distance  the  perform- 
ance of  all  the  arrays  diminishes  to  0 db  since  the  array,  regardless  of 
the  number  of  sensors,  is  in  the  limit  effectively  reduced  to  a single 
probe  (i.e.,  no  resolution  capability).  As  is  apparent  in  Figure  VII-2, 
the  rate  of  increase  in  performance  per  sensor  dwindles  for  larger 
normalized  wavelengths. 

Instrument  noise,  of  course,  degrades  performance  particularly  at 
the  longer  wavelengths.  From  the  dashed  lines  on  Figure  VII-1  which 
represents  individual  array  performance  with  ten  percent  noise  it  ap- 
pears that  the  relative  degree  of  degradation  is  greatest  for  the  densest 
arrays.  In  other  words,  the  denser  arrays  have  "more  to  lose"  in  terms 
of  resolution  capability  and  it  appears  that  increasing  the  number  of 
noisy  probes  is  no  guarantee  of  greater  resolution. 


' ■» 
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VIII.  SUMMARY  AND  ARRAY  DESIGN 


1 . Summary 

The  directional  or  angular  spreading  of  ocean  wave  energy  demands 
spectral  estimation  techniques  of  sufficient  resolution.  Higher  reso- 
lution can  usually  be  achieved  by  increasing  the  size  of  the  array  of 
wave  sensors  though  usually  at  great  cost.  An  alternative  to  greater 
array  size  is  the  nonlinear  MLM  estimation  algorithm.  In  Chapter  III 
it  was  shown  how  1)  directional  spectra  were  naturally  associated  with 
spatial  Fourier  transforms,  and  2)  the  formulation  of  the  direct  and  MLM 
spectral  were  computed  in  terms  of  the  covariance  matrix.  In  Chapter  IV 
it  was  demonstrated  how  the  natural  or  effective  beam  pattern  related  the 
true  spectrum  to  the  estimated  spectrum  in  the  direct  and  MLM  procedures. 
The  interpretation  of  the  beam  pattern  and  spectral  spreading  plots  were 
discussed  in  detail.  Chapter  V presented  the  details  of  the  synthesis  of 
the  directional  spectrum  and  covariance  matrix  used  in  the  simulations. 

In  Chapter  VI,  the  effectB  of  noise,  wavelength,  number  of  elements  and 
spacing  on  array  performance  and  resolution  were  analyzed  using  the 
simulation  of  Chapter  V.  In  particular: 

i)  Noise  contamination  of  the  signal  processing  has  adverse 
effects  on  the  resolution  capabilities  of  the  MLM  technique 
for  all  array  configurations, 

ii)  The  resolution  potential  of  sparse  arrays  increased  dramatic- 
ally with  the  number  of  elements  but  is  not  highly 
dependent  on  the  element  spacing  within  the  array  itself. 

In  Chapter  VII,  a figure  of  merit  was  derived  and  applied  to  the  example 
considered  earlier.  A suitable  figure  of  merit  (FOM)  is  instrumental  in 
the  selection  of  an  appropriate  array  configuration  from  competing 
designs . 

2 . Array  Design 

It  Is  not  yet  possible  to  construct  design  charts  for  the  selection 
of  array  sizes  and  geometry  on  the  basis  of  the  analysis  presented  above. 
A major  unresolved  question  is  the  amount  of  noise  present  in  the  signal. 
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Kst Imntes  of  the  instrument  noise  can  be  made  once  the  sensor  type  has 
been  selected.  Estimates  of  the  amount  of  ambient  field  noise,  both 
propagating  and  non-propagating,  are  more  difficult  to  obtain.  It  will 
therefore  be  necessary  to  assume  conservative  values  for  the  amount  of 
noise  modeled  in  the  simulation. 

It  was  shown  in  Chapter  VI  that  the  MLM  spectral  estimate  approaches 
the  direct  estimate  as  the  noise  level  becomes  large.  This  allows  one  to 
use  the  results  of  the  natural  beam  pattern  analysis  in  order  to  determine 
a lower  bound  on  the  resolution  capabilities  of  a given  array  configuration. 
If  such  a worst  case  design  is  geometrically  and  economically  feasible, 
then  the  actual  MLM  resolution  will  presumably  be  considerably  better.  If 
the  conservative  approach  does  not  yeld  a practical  array,  more  accurate 
noise  estimates  and  additional  simulation  should  be  attempted. 

A simplified  approach  to  the  preliminary  design  or  selection  of  an 
array  for  a specific  application  is  offered  below: 

1)  Estimate  the  range  of  wavelengths  of  Interest.  A sample  point 
spectrum  will  indicate  the  frequency  range  containing 
significant  ambient  wave  energy.  The  particular  ap- 
plication of  the  data  will  generally  suggest  a maximum 

and  minimum  wavelength  bounding  the  range  of  major 
Interest.  For  example,  if  the  problem  is  to  predict 
the  motions  of  a ahip  in  a seaway,  the  high  frequency 
limit  might  be  derived  from  a consideration  of  the  roll 
response.  Waves  on  the  order  of  the  ship's  beam  or 
smaller,  for  example,  will  have  little  effect.  The 
low  frequency  limit  should  definitely  be  below  the  region 
of  maximum  motion  response  and  hopefully  include  the 
most  energetic  portions  of  the  wave  spectrum. 

2)  The  next  step  is  the  estimate  of  the  amount  of  resolution 
considered  to  be  essential.  A convenient  approach  is  to 
select  a maximum  acceptable  width  of  the  main  lobe  of 

the  beam  pattern,  2(8j/2).  This  is  indicative  of  the 
amount  the  true  spectrum  will  be  smeared  out  in  the  pro- 
cess of  estimation.  An  upper  limit  of  the  RMS  array  radius 
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can  be  obtained  from  the  natural  beam  pattern  analysis, 
Eq.  (4.8)  as 

18.2  “X 

2MS  C 0 

where  i / 2 ^ half-power  beam  width  divided 

by  two.  The  MLM  beam  width  is  certain  to  be  narrower, 
hence  better  resolving,  than  the  direct  analysis. 

3)  The  short  wavelength  limit  can  be  used  to  establish 
an  upper  limit  for  the  gage  pair  distances.  In 
order  to  avoid  serious  spatial  aliasing,  the  majority 
of  the  rjj's  should  be  less  than  the  minimum  wave- 
length for  which  a directional  spreading  estimate  is 
desired . 

4)  The  most  important  question  concerns  the  number  of 
elements  or  sensors  that  will  be  necessary.  The 
relative  improvement  in  resolution  for  an  Increase 
in  N is  shown  in  Figure  VI 1-2.  One  should  bear  in 
mind  that  the  rate  of  improvement  in  resolution  de- 
creases as  the  amount  of  sensor  noise  increases. 
Probably  the  best  approach  would  be  to  conduct  a 
sensitivity  analysis  for  various  amounts  of 
spreading  in  the  spectral  estimate.  This, 
combined  with  a reasonable  noise  estimate  and  the 
MLM  simulation,  should  provide  a sound  estimate  for 
the  minimum  number  of  elements. 

5)  Array  shape  is  likely  to  be  somewhat  less  important 

than  the  above  estimates  of  rDUC,  r . , , and  N. 

Knb  1 JH13X 

Local  structure  to  support  the  array  and  questions 
of  safety  of  the  sensors  usually  dictate  to  some 
extent  the  range  of  geometries  available.  The 
above  analysis  has  always  assumed  that  the  waves 
could  come  from  any  direction  or  that  the  orientation 
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of  the  array  could  not  be  controlled.  If  this  is 

not  the  case,  array  geometries  that  are  not  essentially 

equally  distributed  are  conceivable. 

In  conclusion,  the  MLM  offers  greater  directional  resolution  poten- 
tial than  any  of  the  direct  methods  of  spectral  analysis.  Signal  noise 
is  inherent  in  all  measuring  situations.  The  quality  of  the  MLM  spectral 
estimate  can  be  seriously  impaired  by  such  noise  contamination.  There- 
fore, it  is  essential  that  accurate  estimates  of  the  intrinsic  noise 
levels  in  the  measuring  system  be  obtained  in  order  to  insure  adequate 
resolution  with  the  smallest  possible  array. 
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APPENDIX  A 
Array  Geometries 

This  Appendix  includes  representative  small  array  geometries. 

Figure  A- 1 shows  the  array  configurations  to  the  same  scale.  The  angle 
6 is  measurad  counter-clockwise  from  the  x-axis. 

Table  A.l  lists  the  coordinates  (ft.)  of  the  array  elements  used  in 
the  simulation  analysis.  Table  A. 2 gives  a number  of  array  geometry 
statistics  for  the  arrays  of  Table  A.l. 


TABLE  A. 2 
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APPENDIX  B 

Directional  Spreading  Models 


The  following  spreading  functions  have  been  used  as  examples  in  the 
simulation  described  in  Chapter  V,  Section  1. 


Omnidirectional  Spectrum 

This  spectrum  represents  a ring  of  random  plane  waves  propagating 
uniformly  in  all  directions.  The  spreading  function  is  flat  over  the  2 T 
radians  of  azimuth,  i.e.,  F(8)  = 1/2TT,  and  the  Fourier  coefficients  from 
Eq.  (5.6)  are: 


1/2. 7T 

O 


n - o 
n 


(B.l) 


The  covariance  function  becomes  simply, 

CCcoJCN)  = S(oS)  JeCWr')  / Ztt  <b-2> 


This  function  may  be  used  to  model  the  phenomenon  of  propagating  noise 
as  well  as  ocean  spectra. 


Raised  Cosine  Directional  Spectrum 

This  spectral  spreading  function  is  given  by 


F(B') - c • cos^° ( 6- Go')  \ 9 - 9c\  < ^ (B>3) 

where  p is  arbitrary  and  "c"  is  a normalization  constant  such  that 
fl(i  * 1 / 2 7T  . This  spectrum  represents  energy  propagating  to  a point  from 
a half-space  about  the  ray  8 ■ 8^.  The  higher  the  value  of  p the  more 
heavily  concentrated  is  the  energy  about  8 - 0^.  It's  expansion  coef- 


-77- 


a. 


je  cosCnB')  cosf3  CG^J© 

ZTf  ^ cod^CS^dQ 


, (ruo^l^.-V 


(B.4) 


A typical  example  of  the  raised  cosine  spectrum  is  the  "cosine  squared" 

spreading  function  (p=2  in  Eq.  B.4)  centered  about  9=0  =0.  In  this 

o 

special  case, 

- A ‘Sin  ( n / 2 ') 

, (t\  = Oj±\  dbZ...,y 

(4-  n1')  ' J ’ ■’  (B.5) 

In  the  digital  computations,  it  was  found  expedient  to  integrate 
Eq.  (5.3)  numerically  for  each  p of  interest  rather  than  compute  the  a 
coefficients  and  then  evaluate  all  the  Bessel  functions  as  dictated  by 
Eq.  (5.7). 


Delta  Function  Directional  Spectrum 

This  spectrum  represents  energy  propagating  in  a single  direction, 
say  along  all  lines  parallel  to  the  ray  0 = 0°.  Mathematically,  the 
spreading  function  is  represented  by: 


= S(e-6^ 


(B . 6) 


and  substituting  into  Eq . (5.3)  yields  directly 

\ brcos(0o~<^ 

C<“>&  “ S * * 


(B . 7) 


This  is  easily  extended  to  an  arbitrary  number  of  spikes,  i.e., 

1 * 

m»i  (b.8) 

M 

where ^ It  fellows  that  the  covariance  function  for  a set  of  M 
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